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We find exact, analytic solutions of the Dirac equation for a charged, massless fermion in the 
background of a charged, dilatonic black hole in AdSs. The black hole descends from type IIB 
supergravity, where it describes D3-branes with equal angular momenta in two of the three inde- 
pendent planes of rotation orthogonal to the world-volume. The Green's function near the Fermi 
surface for a strongly coupled fermionic system can be extracted holographically from an exact so- 
lution of the Dirac equation at zero frequency but nonzero momentum. There can be several Fermi 
momenta, and they take the form &f — q — n — 1/2 (in units of the chemical potential), where q is 
the charge of the spinor, and n is a non-negative integer that labels the Fermi surfaces. Much as for 
holographic Fermi surfaces based on the Reissner-Nordstrom-AdSs solution, the dispersion relation 
of the excitations near the Fermi surface is determined by the geometry close to the horizon, and 
one can obtain Fermi liquid, marginal Fermi liquid, and non-Fermi liquid behaviors depending on 
the value of fcf. 
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I. INTRODUCTION 

Charged black holes in asymptotically anti-de Sittter 
(AdS) spacetime can be regarded as the gravitational 
dual description of certain strongly interacting fermionic 
systems at finite charge density, such as non-Fermi liq- 
uids pp. This is an application of the gauge/gravity du- 
ality which allows us to calculate fermionic Green's 
function by solving the bulk Dirac equation [5HZ] . A par- 
ticularly well-studied example is the Reissner-Nordstrom 
(RN) black hole in AdS as the geometry [8HTT]. We will 
focus instead on a particular dilatonic black hole in AdSs, 
explored originally in Ref. [IS]; see also Ref. [13] for re- 
lated work. 

The dilatonic black hole in question is sometimes re- 
ferred to as the two-charge black hole. From a five- 
dimensional point of view, this is because two of the three 
mutually commuting U{\) subgroups of the SO (6) gauge 
group of maximal gauged supergravity are nonzero and 
equal, while the third is zero. From a tcn-dimcnsional 
point of view, this black hole describes N coincident D3- 
branes with equal, nonzero angular momentum in two of 
the three independent planes of rotation orthogonal to 
the D3-brane world-volume. The dilatonic black hole en- 
joys several advantages over the better studied RN-AdSs 
black hole: 

• The entropy and specific heat of the dilatonic black 
hole are proportional to temperature at low tem- 
perature, as compared to a nonzero, 0(N 2 ) entropy 
at extremality for the RN-AdSs black hole. 

• Exact information about the position and proper- 
ties of Fermi surfaces is available for the dilatonic 
black hole, for massless bulk fermion actions with 
no Pauli couplings. This stands in contrast with 
the RN-AdSs black hole, where one must resort to 
numerics to find kp. (This is even true of Ref. [14], 
in which numerical work led to strong evidence that 



the Fermi momenta are simple algebraic numbers.) 

• Pair creation of fermions near the horizon, and 
back-reaction of the resulting fermionic matter, 
must distort the RN-AdS 5 geometry to some ex- 
tent. But for the dilatonic black holes there is some 
evidence, to be explained below, that pair creation 
of fermions is suppressed. 

A notable disadvantage of the dilatonic black hole is that 
its extremal limit — which will be our main focus — has a 
naked singularity. Any nonzero temperature cloaks the 
naked singularity with a horizon, but as temperature is 
taken to zero, the dilaton as well as curvature invariants 
become larger and larger at the horizon, until at zero tem- 
perature they diverge. Nevertheless it is straightforward 
to pick out physically reasonable boundary conditions for 
fermions: In particular, for oj = one can simply demand 
that the allowed solutions are regular as the naked sin- 
gularity is approached. 

The main aim of the current work is to solve the mass- 
less Dirac equation, 

y»(V p - <gA M )* = 0, (1) 

in the extremal limit of the dilatonic black hole back- 
ground, and to show that the corresponding Green's func- 
tion exhibits one or more Fermi surfaces if q > 1/2. For 
1/2 < q < 1, there is only a single Fermi surface, and Vf 
is not well-defined. For 1 < q < 3/2, there is still only a 
single Fermi surface, but vp is well defined. For q > 3/2, 
there are additional Fermi surfaces at fcp = q — n — 1/2, 
where n is a positive integer. The outermost Fermi sur- 
face has the simplest properties: assuming q > 1, the 
Green's function near the Fermi surface takes the form 



— u> + vp{k — kp) — S(w, kp) ' 
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Formulas generalizing Eq. ([3| to Fermi surfaces with n > 
can be found in Sec. IIVI 

The organization of the rest of this paper is as follows. 
In Sec. [nj we solve the Dirac equation at u = in terms 
of hypergeometric functions, and find the normal modes 
that determine the location of Fermi surfaces. In Sec. |III| 
we study the near horizon geometry (hereafter IR for 
infrared), solving the Dirac equation and obtaining the 
IR Green's function. In Sec. |IV[ we obtain the Green's 
function near Fermi surface by matching the IR solution 
to a zero-frequency solution away from the IR. In Sec.|Vj 
we numerically solve the Green's function and explain 
the main features at general uj. In Sec. |VI| we conclude 
with some discussions. 



II. NORMAL MODES 

The two-charge black hole in AdS$ is determined by 
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which is from a consistent truncation of the type IIB 
supergravity with three U(l) charges Qi = Q2 = Q and 
Q 3 = 0. The solution in the extremal case is 
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The "horizon" for this black hole is at r = 0, which is a 
spacetime singularity. For the non-extremal case, and its 
ten-dimensional lift, see Ref. [T2"] . 

We will solve the Dirac equation for a massless spinor 
in the above background, but we keep the mass term at 
first. If the metric is diagonal and depends only on the 
radial coordinate r, the Dirac equation can be simpli- 
fied by using the rescaled spinor 

* = (-.9.9 rr ) 1/4 *- The 



equation of motion for \1/ is 

[ 7 "(0 M - iqA^) - m]tf = 0. 



(6) 



We assume that the momentum is in the x direction. 
By plugging a single Fourier mode \P ~ e -iut+tkx-^ ^ Q 
Eq. (ml), the equation for is 



+ iy/g^'f-k - m] # = 0. (7) 
We choose the following gamma matrices for AdSs: 
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Then Eq. ^ reduces to two decoupled equations 



+ {-l) a y/^icr 2 k-rn\iP a =0, (9) 

where ipi and tp 2 are two-component spinors. The equa- 
tion for ip 2 is related to the equation for ipi by k — k. 

The (massive or massless) Dirac spinor ^> in the AdSs 
maps to a chiral spinorial operator at the bound- 
ary [SHZ]- The asymptotic behavior of i> a near the AdS 
boundary is 



ipa 



-> a n r 



b a r~ 



(10) 



The expectation value of the boundary spinorial oper- 
ator dual to the bulk spinor "3/ has the form (0*) = 
(0, &i, 0, b 2 ) T . In fact, Cj, = |(1 — 7-)©^, which means 
that the boundary spinorial operator is left-handed. By 
imposing the in-falling boundary condition at the hori- 
zon, we can obtain the retarded Green's function as 
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G a — 



b a 



(11) 
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Note that if we use the alternative quantization, the 
Green's function is G a = —a a /b a , and the boundary 
spinorial operator is right-handed. If m = 0, G\ and G 2 
are related by G 2 = — 1/Gi [TU]; therefore, the alterna- 
tive quantization for G\ is the standard quantization for 
G 2l and vice versa. By taking into account both G\ and 
G 2 , the alternative quantization gives the same Fermi 
momenta as the standard quantization does, when m = 

ma. 

We will focus on ipi = (ui,u 2 ) T in the following. The 
square roots in the Dirac equation can be eliminated, 



following a method which has appeared, for example, in 
Ref. [15J . Define u± = u\ ±iti2. From Eq. ([9]), we obtain 
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A(r) = *W J ^ L (w + gi4 t ), 

(14) 

The Eqs. (12 1 and (13) can be decoupled to obtain two 



second-order differential equations: 



u + + p( r W + + q(r)u + = 
u"_ + p(r)u'_ + q(r)u- = 0, 



(15) 
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After we solve Eq. ( |16| _for u- , we need to plug in u_ to 
Eq. (13) to obtain u+r For the metric we consider, 
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In the following, we will only study the m — case, in 
which p(r) and q(r) are rational functions of r. We are 
most interested in the following two questions: whether 
there are Fermi surfaces, and whether there are quasipar- 
ticles near the Fermi surfaces. We will solve the Dirac 
equation at u — first, and the solution indicates that 
there are one or more Fermi surfaces when q > 1/2, as 
summarized more precisely in the text following Eq. (JlJ . 
Then the perturbation at small lu will give the Green's 
function near the Fermi surfaces. 

When uj = 0, the boundary condition for at the horizon 
is that the solution is regular. The solution for u± can 
be written a£] 
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FIG. 1. The real axis in the complex r-plane maps to a 
circle in the complex z-plane. The hypergeometric function 
F(a, /3; 7; z) has a branch cut from 2 = 1 to 00 in general, but 
the branch cut is absent when a is a non-negative integer. 



where 



V2Q' 



(21) 



The chemical potential V%Q is a unit of the energy scale. 
To have physical bound states, Q and q must have the 
same sign; we assume Q > and q > 0. This system has 
rotational invariance; we can choose k = (A;, 0,0), where 
k > 0. Thus we have v k > 0, without loss of generality. 

By defining v k — q+ 1/2 = — n, the solution for u\ and 
u 2 is 



Ul = 



u 2 



(-1)* 



(-1)' 



2i 2i 
The Green's function G(w, k) at to = is real 

,(-l) n+1 u*_ - u_ 



(22) 
(23) 



G\ = lim — 

r->oo 



= lim — 



{-\)n+l u *_ +u _ 



= G\. 



( 24 ) 

This apparently implies that the spectral density is zero 
at oj = 0. However, we need to shift the pole at u> = 
by w — > w + ie, and then we will obtain a delta function 
in the imaginary part[^] 

The normal modes are determined by Uil^oo = 0. In 
general, the hypergeometric function 2-F1 {ot, /?; 7; z) has 
a branch cut from z = 1 to 00. At the AdS boundary, 



2 F 1 (-n, v k ; 2v k + 1; 2 - ie). 



(25) 



We cannot only solve Eqs. ( |15| and | |16[ l and discard Eqs. \12\ 

and |TJ|, because there are 6riy two bWndary conditions. After For example, for a free electron near k F (k ± = k - k F ): 

we solve u- from the second-order equation is fully 1 1 . , . 

determined by the first-order equation \13) . — w vpkj_ — ie ~ —u> + vpkj_ ^ VF 

Note that (-1)" := (-1 + = e ina and (-1 - ie) a = e~ iva . 



4 



Thus, U- and u*_ take values at different sides of the 
branch cut, as shown in Fig. [T] However, if a = —n, 
where n = 0, 1, 2, • • • , the hypergeometric function is 
an nth-order polynomial of z, and the branch cut from 
z = 1 to oo is absent. More generally the equation that 
determines the normal modes is Eq. (A7| in appendix [Aj 
in which we conclude that there are no physical solutions 
when n is not a non-negative integer. 

At the AdS boundary r — > oo, u*_ = u_ if a = —n, 



where n = 0, 1, 2, 



Therefore, if n is even, 



An) _ 



q — n — 1/2 gives the Fermi surface for the standard quan- 
tization (ui = 0); if n is odd, = q — n— 1/2 gives the 
Fermi surface for the alternative quantization (u 2 = 0). 
This conclusion is for the Green's function Gi, which ob- 
tained by the upper-half components of the bulk spinor. 
Recall that G 2 — — 1/Gi, which is obtained by the lower- 
half components of the bulk spinor. In the following, we 
use the standard quantization only. Taking into account 
both G\ and G 2 , we conclude that the Fermi momenta 
are determined by = q — n — 1/2, where n is a non- 
negative integer such that q — n— 1/2 > 0. Note that the 
alternative quantization gives the same Fermi momenta, 
with the difference that the boundary fermionic operator 
is right-handed. 

By perturbation, we can obtain the analytic solution of 
the Green's function near the Fermi surface. The Green's 
function can be written as 



G R (u,k) 



-uj + vp{k — kp) — S(w, kp) ' 



(26) 



where kp = |kp| is the Fermi momentum, vp is the Fermi 
velocity, and 



k) = h{k)Q k (uj), Gk{u) = c(k)uj 



2"k 



(27) 



As fermionic Green's functions, G and Q satisfy Im(G) > 
and Im((?) > for all real u). The result shows that 
vp > 0, Z > 0, and h > 0. The Fermi momenta are 
determined by 
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(28) 



where n = 0, 1, 2, • • • , \q— 1/2J . When n is even, Eq. (26) 



is for Gi; when n is odd, Eq. (26) is for G 2 . Again note 
that we need to shift the uj — pole to the lower half 
complex w-plane by uj — > to + it to obtain a well-defined 
retarded Green's function. 
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FIG. 2. The inner (near horizon) and outer regions, where 
the solutions of the Dirac equation are denoted by ipi and ipo , 
respectively. 



treated as small perturbations no matter how small uj is. 
This section and the next are in parallel with Ref . |10j , in 
which a systematic method is developed for treating the 
extremal black hole system. Usually this method relies 
on numerics to fix certain quantities, such as the Fermi 
velocity. The example we provide is exactly solvable, in 
the sense that a perturbative treatment of the small uj 
regime can be obtained through matched asymptotic ex- 
pansions of analytically known functions. 

We divide the geometry into inner and outer regions, 
as shown in Fig.[2j The inner region refers to the IR (near 
horizon) geometry, in which the Dirac equation can be 
exactly solved to give an IR Green's function. The outer 
region refers to the remaining geometry, in which we can 
make perturbations for small uj. Then we need to match 
the inner and outer regions. 

The IR geometry is examined as follows. In the r — » 
limit, the metric becomes 



ds 2 
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Therefore, the IR geometry is conformal to AdS 2 x M 3 . 
This can be made more explicit by change of variables 
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and the metric becomes 
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The gauge field A t becomes 



$ = 



QC 2 



(30) 



(31) 



(32) 



III. IR GEOMETRY AND GREEN'S FUNCTION 

We expect that the Green's function near the Fermi 
surface can be obtained by the perturbation of small uj 
around the exact solution. However, the ordinary per- 
turbation method is not enough when the black hole is 
extremal. As pointed out in Ref. |10) . when it is suffi- 
ciently close to the horizon, w-dependent terms cannot be 



We will switch back to the r coordinate. Note that in 
the RN-AdS black hole system, $ ~ r, and thus the 
electric field E = V<I> is constant at the horizon. In our 



system, $ 



and thus the electric field E = V<I> 



falls off toward the horizon. This leads to a significant 
difference relative to the Dirac equation in AdS 2 . In the 
near horizon limit r — > (£ — > 00), the contribution 
by the electric field to the Dirac equation is negligible. 
Nevertheless, the flux is conserved by d(e 4a *F) = 0. 
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We solve the Dirac equation in the geometry Eq. ( 29 ) 



without the electric held. The solution for u± with in- 
falling boundary condition^] is 

u- =C^W 1/2> „ h {-iu/r) (33) 
u + =ii/Cy/fW- 1/2tVh (-iu}/r), (34) 



the black hole horizon [TUl |T7j . It has been argued that 
back-reaction from the pair production alters the IR re- 
gion to a Lifshitz geometry [18 ; a candidate of the hnal 
geometry was constructed as the electron star [T5rt21| . 
In our system, v k is always real, and the electric held 
approaches zero in the near horizon limit. 



where W is a Whittaker function, and C is a constant. 
Denote ipj as the solution in the inner region. In the 
near boundary limit of the IR geometry, the asymptotic 
behavior is 



ipi — > a 



U!\ v k 



as 



r 



0. (35) 



More precisely, the inner region solution can be written 
as 

^ = i)/ l (l + -)+ftM»-f"" l (l + •••). ( 36 ) 

where v± must be chosen to match the normalization of 
Eqs. (19) and (|20[ ). By expanding t/jj — (111,1*2) from 
Eqs. (33 1 and 4346, we know that v± take the following 
form 



Aj 



A_ 
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(37) 



where A± are constants. The Green's function Q k {ui) de- 
pends on the ratio A + /A_. However, the self-energy E 
is independent of A± after matching the inner and outer 
regions. We choose A + = A_, and then the IR Green's 
function i^£] 



g fcM = e ^(V2-, fc) r(i/2-^) 



(38) 



The IR Green's function can be generalized to hnite tem- 
perature when the back hole is near extremal: 
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The main difference between the RN-AdS black hole 
system and our system is attributed to the IR geometry 
with the gauge held. In the RN-AdSs black hole system, 
the IR geometry is AdS2 x IR 3 , and the electric held is 
nonzero at the horizon . The IR scaling exponent has the 
form v k = \Jh 2 — k%, which depends on the charge of 
the spinor, and will become imaginary if the charge is 
large. The system with this IR behavior is studied as 
a semi- local quantum liquid |16j . The imaginary v k im- 
plies an instability causing by the pair production near 



IV. GREEN'S FUNCTION NEAR THE FERMI 
SURFACE 

In the outer region, the solution at small u can be 
written as 



where 



^0 = v+ + Gk(u)y-, 



(0) , (1) , 2 (2) 

r}± = r]± + uir]± +uj r)± 



The asymptotic behavior near the horizon is 

r -> 0, 



rj® = v ± r ±v " 



(40) 
(41) 
(42) 



which is matched with the inner region solution, Eq. ([36| . 

Here 77^°' = (ui,U2) T , where u\ and 112 are solutions in 
the outer region as Eqs. (22 1 and (23). We expand Ui in 
the r — > limit 
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-r Uk (! + ■■■)■ 



(43) 



Similarly, we can expand ui and the solution of rji^ . The 
normalization constants v± are 



v± 



±1 



The asymptotic behavior near the boundary is 



(44) 
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Consequently, the Green's function near w = to the 
hrst order is [TU] 
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We only summarize the result of the perturbation 
method given in Appendix C of Ref. |10) . Some nota- 
tions are slightly changed here. Dehne 



4 The in-falling wave in terms of the coordinate f is e 1 "^ as £ — > 00. 

5 Other ways to write down Gfe(^) are 
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J x = (*o,r a! *o) 



drVft^^fjVa^f, (47) 



which are integrations of hypergeometric functions in our 
system. Note that both J* and J x are negative. The 
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various functions in the Green's function Eq. (26) are 
determined as follows: 
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where the quantities above are evaluated at k — kp. 



(48) 



By plugging u\ and into Eq. (47), the integration 



can be evaluated for non-negative integers n. The first 
three results of J* are 
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By induction, we find that the nth J* is given by 
jt (n) = 

n!0r]8^ 2 + {An + 2)v - l]T{v + l/2)Y{v - 1/2) 

2™+ 3 v / 2Q(2i/ + l)r(iv + n/2 + l/2)T{v + n/2 + 1) ' 

(50) 

where v = vY 1 ' . 

After we obtain J x , the Fermi velocity is 



(n) 2(2»/ + l)(2i/-l) 
' F 8^ 2 + (4n + 2)^- 1' 
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where ^ = . When n is even, is for G?xi when 

n is odd, is for G2. If we take V2Q = 1, the only 
independent parameter is the charge of the spinor, in 
terms of which the vp can be written as 



(n) = 4(g- n)(q- n- 1) 
Jp 4q 2 - 3(2n + l)q + 2n(n + 1) ' 



(52) 



We can see that < Vp < 1, and vf — > 1 as q — > 00. The 
Fermi velocities as a function of the charge is plotted in 
Fig-! 

If n is even, for the Green's function Gi, 



, (n) _ 2^2Qr(n/2 + l/2)r(i/ + n/2 + 1) (rl) _ 
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if n is odd, for the Green's function G2, 
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(53) 



(54) 



where v 
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The ratio Z/vp as a function of charge 



is plotted in Fig. [4] The self-energy is given by 

(n) = r(2^ + n + l)r(l/2 - i/)e fa (V 2 -^a^ (TO) 
2 6l '- 1 (y2Q)^- 1 r(n + l)r(^ + l/2) 3 ' 
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FIG. 3. Fermi velocity as a function of charge q, where 
q > n + 1 for the nth Fermi surface. The solid lines are for 
n — 0, 2, • • • , and the dashed lines are for n = 1, 3, • • • (from 
left to right). 




FIG. 4. Z/vp as a function of charge q, where q > n + 1/2 
for the nth Fermi surface. When q — n + 1/2, Z/vf = 2/7T, as 
indicated by the horizonal line. The solid lines are for n — 0, 
2, • • • , and the dashed lines are for n — 1, 3, ■ • ■ (from left to 
right). 



where v = 



The spectral density p = Im(G) as a 



function of ui at different values of k is plotted in Fig. [5] 
As k — > /cp, the quasiparticle peak will become a delta 
function at k = kp. 

Another way to write down the Green's function is 



G{u,k) 



hi 



(56) 



k± u) 

where k± = k — kp, hi = Z/vp, hi = \Yi / (<jj 2v v p)\, and 
7fc = tt(1/2 - v k ) + ar g r(l/2 - v h ), (57) 



The Green's function in the form of Eq. ( 56 1 is analyzed 



in Ref. [TU] in detail. The poles never appear in the 
upper half complex cj-plane of the physical sheet. The 
three cases v k > 1/2, — 1/2, and z>% < 1/2 correspond 
to Fermi liquid, marginal Fermi liquid, and non-Fermi 
liquid, respectively. 

The Green's function for the non-Fermi liquid {y% < 
1/2) can be written as 



G = - 



Cjk 1 / 2 ^- 1 
u-czk 1 / 2 "* 



(58) 
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FIG. 5. Spectral density as a function of ui at three different 
values of k. We plot the near u = region for q = 4.5 and 
n — 2, which give fcf = 2. The dotted, dashed, and solid 
curves are for k = 1.8, 1.9, and 2, respectively. When k — Jcf, 
the quasiparticle peak becomes a pole. 



FIG. 6. Spectral density as a function of u) at k = 2 with q = 
10 and \/2Q = 1. We can see five peaks for the quasibound 
states. As we increase k from small k, the peaks will move 
to the right. Each time a pole go across ui = 0, we obtain 
a Fermi momentum. Therefore, there are five Fermi surfaces 
from Gi , as the formula of k^ predicts. 



where 



hi 



2v{h 2 e i ~ ik ) 1 l' 2vk 



(>2 



(/i 2 e i T fc ) 1 /2^ 



(59) 



The residue vanishes as k — > kp. The pole cj* moves 
along a line approaching the origin with the angle 



arg(w.)= \ )l ' 4 ?\ 



k > kp. 



(60) 



We can see that 0* ^ (0, 7r), which is the upper half 
plane of the physical sheet 6> g (— ir/2, 3tt/2). There is a 
particle-hole symmetry due to lm(uj*)\k ± = lm.(ui*)\-k±- 
What is especially interesting is the marginal Fermi 
liquid. In the v k — > (l/2) + limit, 

= -u - u(2\nuj - in + c^)e + C(e 2 ), (61) 

where e = — 1/2, and is a real constant. We 
can see that the —uj in E exactly cancels the — to in the 



denominator of Eq. ( 26 ) , which is a delicate cancellation 



between the UV and the IR data. The Green's function 
for the marginal Fermi liquid is 



G 



hi 



k±_ + \ (n + l)u\nu + cWcd' 



(62) 



where hi and can be easily obtained by the exact 
solution. The residue also vanishes as k — > kp. 



V. GREEN'S FUNCTION AT ARBITRARY u 

To obtain the Green's function when u is not small, 
we can solve the Dirac equation numerically with the 
boundary condition near the horizon as Eqs. (33) and 



(34 1. Alternatively, we will solve the flow equation for 



£ = u-ijui as follows 



2m 



fip(LU + qA t ) 




(63) 



with the boundary condition £| r= o — i ^ 0). The 
Green's function is obtained by r 2m £| r _ i . 00 . A typical 
case of the spectral density p — Im(G) as a function of 
uj is plotted in Fig. |6j in which the peaks correspond to 
quasibound states. 

The non-analytic features of the Green's function from 
the RN-AdS black hole at finite temperature are studied 
in detail in Ref. [22] ■ The poles of the Green's function 
are schematically plotted in Fig. [7J in which we ignore a 
small difference that the poles in the finite temperature 
case cannot be exactly at the origin u = 0. The RN- 
AdS black hole system and our system have some similar 
features, as follows. Consider the m — case. There are 
no poles in the Green's function at k = 0, so we start 
with a small k. If the charge of the spinor q is sufficiently 
large, there are quasibound states. As we increase k, the 
poles with Re(w) < will move to the right. When a 
pole goes through the origin u) = 0, we obtain a normal 
mode, which indicates a Fermi surface. The number of 
the quasibound states equals the number of the Fermi 
surfaces. 

The schematic plots of the dispersion relation and the 
Fermi surface are shown in Fig. [5J Our system has at 
least one Fermi surface when q > 1/2. If we increase q, 
more Fermi surfaces will appear, and the Fermi surfaces 
are equally spaced. 

The number of quasibound states can be estimated by 
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FIG. 7. Schematic plot of the poles of the Green's func- 
tion. The generic feature is that there are quasibound states 
when the charge of the spinor is sufficiently large. The highly 
damped modes are plotted from the RN-AdS black hole sys- 
tem. 




FIG. 8. Schematic plot of dispersion relation and the Fermi 
surface by the massless spinor in the bulk (for massive spinor 
in the bulk, see Ref. [22]). The system has rotational in- 
variance, and we only show some intersections. If we decrease 
the charge q, the Fermi surfaces will shrink but keep the same 
space. 



the WKB method. The effective potential i^] 



v:, 



9 —k 2 

g rr 



9' 



(64) 



The distinctive shapes of the effective potential are plot- 
ted in Figs. [9] and [To] At the AdS boundary, V cS = 0. 
The near horizon behavior is 



uj 



k 2 



quj 

7^Q 



UJ 

W 

m 2 



2 Q2/3 r 4/3 



(65) 



We need to treat uj = and uj ^ cases separately. 
Assume at least one of k and m is nonzero, otherwise the 
Green's function has no poles. 

When u = 0, the leading term in V e g is positive as 
r — > 0, which implies that the state cannot tunnel to the 
horizon and thus is stable. This is similar to the electron 
star, but different from the extremal RN-AdS black hole, 
in which there are no stable bound states when Vf. is 
imaginary. 

When uj 0, the leading term in V Q g is negative, which 
implies that there are no exact bound states with uj =/= 0. 
The state can tunnel through a barrier to the horizon, 
which will lead to an imaginary part of the modes. The 
qualitative features of the effective potential are similar 
to the RN-AdS black hole. If q is large enough, there is 
a potential well with a barrier. The quasibound states 
in the well can tunnel through the barrier. In Fig. [7j the 
modes near the real uj axis correspond to the quasibound 
states. 



VI. DISCUSSION 

Starting from a dilatonic black hole derived from a 
consistent truncation of type IIB supergravity, we have 
studied the fermionic Green's function dual to massless 
fermions in the bulk. We obtained exact analytic results 
at zero frequency, and exact asymptotic results for small 
frequencies. These analytic results capture key features 
of the strongly coupled fermionic system modeled by the 
gauge/gravity duality, and they provide a new universal- 
ity class for the strange metal phase at quantum criti- 
cality. Provided that the charge of the bulk fermion is 
not too small, there are Fermi surfaces. Their Fermi mo- 
menta are equally spaced, and there are a finite number 
of them, approximately proportional to the charge of the 
bulk fermion. The IR scaling dimension is always real. 
The properties of the RN-AdS black hole system and our 
system are compared and summar ized in Table [I] 

There are several instabilities that can modify the 
bosonic background, including the superconducting and 
the Gregory-Laflamme instabilities; however, these insta- 
bilities all involve extra fields not present in our consis- 
tent truncation of the supergravity Lagrangian, Eq. Q. 
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FIG. 9. Effective potential when uj = 0. 
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FIG. 10. Effective potential when u) ^ 0. 



TABLE I. Comparison between the RN black hole system and our system, with the same UV geometry as AdSs. 



RN black hole 



two-charge black hole 



Charge 
Entropy 

IR (near horizon) geometry 
Electric field near horizon 
Stability near horizon 
IR scaling exponent Vk 
v k | 1/2 



Qi = Q 2 = Q 3 = Q 
S — constant 
AdS 2 x R 3 
E = constant 

unstable due to pair production 
Vk oc y/k 2 — k 2 



Qi = Qi = Q, Q 3 = 
SocT^O 

conformal to AdS 2 x R 3 
E oc r -> 

stable against pair production 

v k oc k 



Fermi liquid, marginal Fermi liquid, non-Fermi liquid 



Appendix A: Mathematical notes 

Hypergeometric function 2-Fi(a! ) /3; 7; z). We will de- 
note 2-F1 by F for simplicity in the following. The deriva- 
tive of the hypergeometric function can also be expressed 
by a hypergeometric function: 

~F(a, P; 7; z) = °^F(a + 1, /3 + 1; 7 + l;z). (Al) 
dz 7 

The following formula can be used to combine the sum 
of two hypergeometric functions: 

jF(a, P; 7; z) + azF(a + 1,0 + 1; 7 + 1; z) 

= 7 F(a,p + l; T ,z). (A2) 

In general, the hypergeometric function has branch 
points at z = 0, 1, and 00. By convention, we make 
a branch cut from z = 1 to 00, and take the principle 
branch as — 2ir < argz < for \z\ > 1. The follow- 
ing formula can be used to transform a value above the 
branch cut to another value below the branch cut: 

F(a,/3; T , z) = (1 - z)- a F(a, 1 - p; 7; ^-). (A3) 



The z — 2 point has the following special property: 

y 

2 ± ie — > 2 =F ie. 
• " ' T we have 



z - 1 



(A4) 



By Eq. {A3) and (-1 - ie)~ a = e" 

F(a, /?; 7; 2 + ie) = e ia *F(a, 7 - /?; 7; 2 - ie). (A5) 
We define 

F(a,/3; 7 ;2) := ^; 7; 2 - ie). (A6) 

The condition for the normal modes is 

F(-n, v k - 2v k + 1; 2) = ±F(-n, ^ fe + 1; 2^ fe + 1; 2), (A7) 

where n = — v\. + q — 1/2. The plus sign is for G\ and 
the minus sign is for G2. We assume that n is a non- 
negative integer at first. If n is even, the above equation 
with the plus sign is satisfied; If n is odd, the above 
equation with the minus sign is satisfied. We can nu- 
merically check that they are the only solutions when 
q > 0. When q < 0, there is another set of solutions 
due to the q — > —q, uj — > — uj, u\ O U2 symmetry of the 
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Dirac equation; however, these solutions are unphysical 
because they give Im(G) < 0. Intuitively only if a parti- 
cle and the black hole have the same charge can there be 
a balance between the attractive gravitational force and 
the repulsive electromagnetic force on the particle. 
For non-negtive integer n, 



F{-n, v\ 2v + 1; 2) 

f r(n/2 + l/2)r(i/ + l/2) 

y4l> + n/2 + l/2) 
I>/2 + !)!> + 1/2) 

v / 7rr(V + n/2 + 1) 



if n is even 



if n is odd. 



(A8) 



Gamma function. Useful identities for the Gamma 
functions include 



'(» + ;) 



(2n)! (2n-l)H 

7T = — \]-K 



4™n! 
T(z)T(l-z) = 

smirz 

T(z)T(z + 1/2) = 2 1 ~ 2z 0r T(2z) 



(A9) 



Whittaker function. Whittaker's equation is 
d 2 W ( 1 A 1/4 -/i 2 



dz 2 



4 z z 2 



W = 0. (A10) 



We can write down the general solution as CiW\ tll (z) + 
C2W-\^(—z), where for large \z\ one has 

W x Jz)~e- z / 2 z x (l + ---), \z\^oc. (All) 

As special cases, W ± i/ 2 ^{z) are related to the modified 
Bcssel function K v (z) by 

^-v,,w-^(w(i)-w(i))- 

(A12) 
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